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Recent experiments have shown that when a near-hemispherical lipid vesicle attached 
to a solid surface is subjected to a simple shear flow it exhibits a pattern of membrane 
circulation much like a dipole vortex. This is in marked contrast to the toroidal circula- 
tion that would occur in the related problem of a drop of immiscible fluid attached to a 
surface and subjected to shear. This profound difference in flow patterns arises from the 
lateral incompressibility of the membrane, which restricts the observable flows to those 
in which the velocity field in the membrane is two-dimcnsionally divergence free. Here 
we study these circulation patterns within the simplest model of membrane fluid dynam- 
ics. A systematic expansion of the flow field based on Papkovich-Neuber potentials is 
developed for general viscosity ratios between the membrane and the surrounding fluids. 
Comparison with experimental results [C. Vezy, G. Massiera, and A. Viallat, Soft Matter 
3, 844 (2007)] is made, and it is shown how such studies could allow measurements of the 
membrane viscosity. Issues of symmetry-breaking and pattern selection are discussed. 



1. Introduction 

One of the defining features of plant cells is the vacuole, an organelle filled with water 
that acts as a nutrient storehouse for the cytoplasm. The vacuole is contained within the 
vacuolar membrane, a lipid bilayer structure known also as the tonoplast. It is a nearly 
ubiquitous feature of plants that the cytoplasm, the thin fluid layer surrounding the vac- 
uolar membrane, is in constant motion through the phenomenon of cytoplasmic streaming 



( Verchot-Lubicz & Goldstein 2010) in which motor proteins move along filaments and 



entrain fluid. Dating back to important work by Pickard (1972) it has been suggested 



that the shear created in the cytoplasm is fully transmitted through the tonoplast into 
the vacuolar fluid, a conjecture supported indirectly by recent whole-cell measurements 
of the streaming velocity profile in large plant cells ( van de Meent erZ|2010| ). 

In considering the dynamics of shear transmission across the tonoplast, ideas may 
be drawn from the large body of work on the response to flow of lipid vesicles. These 
are closed bilayer membrane structures enclosing fluid, whose membrane is generally in 
the 'fluid' phase with zero shear modulus. These studies constitute a highly active area 
of research complementing the well-established understanding of membranes in thermal 



equilibrium (Seifcrt 19971. The first insights into vesicle motion under flow were made 



when Keller & Skalak (19821 described prototypical 'tank-treading' and 'tumbling' be- 



haviour. Further developments were quickly made in the context of red blood cell dynam- 
ics; Barthes-Biesel & Sgaier ( 1985 ) used a visco-elastic membrane model to understand 



the observed tank-treading behaviour, and Feng et al. (1989) provided detailed insights 
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into the precise structure of the area-preserving membrane flow patterns. However, it is 
more recently that advances have been made in understanding the full three-dimensional 
behaviour of nearly-spherical vesicles ( Seifert||1999 ) , revealing a rich phase space of both 



stable ( |Misbah||2006| |Lebedey et aL|2007[ |Deschamps et a/.||2009| |Zhao fe Shaqfeh|2011 ) 

and unstable (Kantslcr et al. 2007) behaviours. 

Less well-studied has been the hydrodynamics of vesicles when constrained to lie stuck 
to a surface. Such problems arise as the natural extension of considering the flow of these 
objects near to walls, since adhesion can occur upon the vesicle coming into contact with 
the surface ( Abkarian fe Viallat|2 008 ) . Indeed, the elementary problem of an immiscible 
fluid droplet in contact with a no-slip plane immersed in a bulk flow, absent the inter- 



facial membrane of a vesicle, is one that has attracted detailed study. Dussan V. ( 1987 I 



performed a comprehensive theoretical study of the behaviour of shallow droplets in such 
situations, deriving a flow pattern whereby the droplet flows with the applied shear on 
the surface and recirculates back underneath inside. Sugiyama fc Sbragaglia] p008l con- 
sidered a perfectly hemispherical droplet with similar conclusions. The presence of an 
interfacial membrane disallows such flows, since they are not area-conserving and the 
membrane is unable to flow back through the vesicle, and so a qualitatively different flow 
pattern is guaranteed simply by the inclusion of the membrane. 



Lorz et al. ( 2000 ) studied the weak adhesion of vesicles and observed a tank-treading 



type behaviour of the membrane under shear flow, with symmetric recirculating regions 
on either side of the vesicle midline. To further quantify the membrane hydrodynamics, 
|Vezy et al. (2007) studied adhered vesicles in shear flow with varying contact areas, and 
observed the same symmetric doubly-recirculating pattern of motion of the membrane in 
all cases. These experiments have yet to be complemented by any theoretical modelling 
and calculation of the expected membrane and bulk flow patterns in such geometries. 
Therefore, as a first step towards understanding these observations we formulate here a 
simplified theory for a vesicle adhered to a flat substrate subjected to shear flow. The 
vesicle is assumed to be an axisymmetric spherical cap, and the membrane to be an 
impermeable, two-dimensional, incompressible Newtonian fluid, where we approximate 
the bilayer structure of the membrane as a single layer by assuming zero inter-layer slip. 
On the phenomenological length scales of this problem (~ 1 — 100 U-m) the inter-layer 
dissipation is low (Seifert & Langer 1993), and so this simplification is justified. We 



will also assume that the vesicle adheres sufficiently strongly to the surface so that we 
may neglect any shear-driven deformations allowed by remaining excess area. (Typical 
experimental shear rates fall in the range 1 — 10 s _1 ; see |Vezy et aL|[2007| ). Throughout 
we work in the Stokes regime Re <^ 1 for all fluids. 



2. Theory 

2.1. Problem formulation 

The vesicle is a spherical cap of radius R and origin x = y = 0, z = H, adhered to the 
plane z = 0. Let (r, 9, 0) be spherical polars centred on x = y = z = 0, and similarly let 
(s, ip, (f>) be centred on the origin of the adhering vesicle. The vesicle then has a maximum 
extent ip = a at 9 = 7t/2, where a = 71— arccos(£f/i?). Figure [I] illustrates the geometry. 
We solve for three velocity fields: the flow u~ inside the vesicle s < R, the 2D flow u m on 
the membrane s — R, and the flow u + outside the vesicle s > R which we will write as 
u + = U + + jzx, where here and henceforth hatted quantities denote coordinate system 
basis vectors. To these fluids we associate respective viscosities 7/_, rj m , rj + . 

In what follows we non-dimensionalise space by R, velocities by Rj and volume pres- 
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Figure 1: Geometry of the problem. 



sures by j 2 pR 2 . Boundary conditions are then applied on s — 1 and the imposed shear 
rate is unity. We also define h = H/R, the effective height in this scaling. The external 
flows obey the unforced Stokes and incompressibility equations, 

KiVV - Vp* = 0, V-w ± =0, (2.1) 

where K± = r\±j A ipR 2 . We impose far-field asymptotics, planar no-slip and radial no- 
penetration: 



u ~ zx as r — > oo, 
.± 



Oonfl = 7t/2, 
• s = on s = 1. 

All three velocities must be continuous across the membrane, 

,+ _ „,m _ „,— „ ^ 



it 



m on s 



(2.2) 
(2.3) 
(2.4) 

(2.5) 



which, along with (2.3 1, implies the planar no-slip condition u m = at ip = a. 

In the absence of the membrane, we would only impose continuity of the bulk flu- 
ids' normal stresses at the interface, but since membranes can support tension, the bulk 
stresses may be discontinuous. By assumption the membrane satisfies the Stokes equa- 
tions and incompressibility constrained to the surface s = 1. For a planar membrane 
this is simple (e.g. Lubensky & Goldstein 1996), but curvature induces an extra term. 



Henle & Levine (2010) (see also Henle et al. 20081 formulate the Stokes equations in 
terms of covariant derivatives on an arbitrary manifold, from which the membrane ten- 
sion (equivalent to a two-dimensional pressure) may be eliminated by taking the curl, as 
in the familiar derivation of the vorticity equation. In our notation, 

2e+ 

J s— 1 



8 ■ V X 



\V 2 u r ' 



s ■ V x 



2e: 



= 0, V • u m = 0, 



(2.6) 



where = e s ^xjj + e s< p4> are the bulk fluids' in-plane normal rates-of-strain, V denotes 
the gradient operator constrained to the surface s = 1, and we define 

Vm 

r+ — . 

Rv± 

These are the non-dimensional form of the 'Saffman-Delbriick' lengths l± 
( Saffman fc Delbriick|1975 Saffman|1976 Henle fc Levine|2010 |, and are the parameters 
with which we control the membrane dynamics. Using the typical values n± ~ 10 -2 Poise 



(2.7) 
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(water), n r , 

relevant experimental ranges to be r± 



10- 



10. 



2.2. Solution method 

We will first construct a set of basis functions for the bulk fluids such that the planar no- 



slip condition (2.3) is automatically satisfied. A common approach to solving Stokes flow 



problems is to use Lamb's solution (Happel & Brenner 1991), and Ozarkar & Sangani 
( |2008[ ) showed that it is possible to write down a no-slip image system for each individual 
Lamb mode in the manner of the original solution of |Blake| ( |1971[ ) for a simple Stokcslct. 
However, this becomes algebraically unwieldy, and we adopt a slightly different approach. 

A method similar to using a Lamb expansion is to writ e down the solution in terms of 
Papkovich-Neuber potentials ( |Tran-Cong fc Blake||l982| ): if V 2 A = and V 2 P = 0, 



v = V(r • A + B) - 2A, p = 2rjV • A, 



(2.8) 



is a solution to the incompressible Stokes equations r/V 2 v = Vp, V • v = 0. This reduces 



the problem to one of solving the vector and scalar harmonic equations. Shankar ( 2005 ) 



uses this representation to solve for Stokes flow inside a circular cone using three harmonic 
basis functions, one scalar and two vector. The scalar harmonic is the usual solution to 
Laplace's equation in spherical coordinates, 



B(r; i>, m) 



b P™(cos8), 



(2.9) 



where v G M (since the solution need not be analytic in z < 0) and m£Z, and P™ are 
the generalised associated Legendre functions. Two independent vector harmonics can 
then be derived from B as A\ — r x V(rP) and A 2 = Bz, which evaluate to 



Ai(r;v,m) = -r v e im *[m 
A 2 (r; m) = rV m< ^ [cos 0P™ (cos 



csc6P™(cos60<? + iP; 

- sin 0P™ (cos 8)0], 



,m '(cos6»)0], 



(2.10) 
(2.11) 



where primes indicate differentiation with respect to 0. Since the only non-zero imposed 
velocity is simple shear flow we only need the m = 1 mode, and can take the real part of all 
basis functions, whose m-dependence we henceforth omit. We also write V u {x) = P„(x). 

The next step is to find combinations of these harmonics such that the no-slip condition 
(2.3) is satisfied. In Shankar ( |2005 ) this is done by evaluating the velocity field on the 
cone using one unit of B(r; v + 1), a units of Ai(r; v) and b units of A 2 (r; v). However, 
the case of a flat plane is a special one and solutions with no contribution from B are 
possible, so here we must take an arbitrary c units of B(r; v + 1); then (2.8) gives 



v r = r v cos (j)[c(v + l)T> v+ i+b(v -1) cos 9V V ], (2.12a) 
v e = r v cos <p [cP' v+1 + 2a esc 6 V u + b {sin 6 V„ + cos dVl}] , (2.126) 
= — v v sin (f) [c esc 9 *Pv+i + 2aV' v + bcotdV,,} . (2.12c) 

(From here on, where unspecified the argument of V v is cos#.) It is simple to find that 
demanding v — at 8 = n/2 leads to two possible conditions: 



Either 2a + b - {v + 2)c = and V u +i (0) = 0, 
or a = c = and V v {0) = 0. 



(2.13) 
(2.14) 



The condition P„ (0) = has solutions v = 0, 2, 4, . . . and v = —1, —3, —5, . . .. (Note that 
V = P_i = 0.) 
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With this basis the general solution with no-slip on z = can be written 

oo 

v(r,6,<f>) = V r u [a„F v (9,cl>) + b v G 1/ {6 : <p) + c u H„(6, </>)], 



(2.15) 



where the functions F,G,H may be read off from (2.121 as the coefficients of a,b,c 
respectively. The coefficients must satisfy 



2a„ 



u -(v + 2)c„ = for v = 1, 3, 5, . . . , v = -2, -4, -6, . . . , 
a v = c v — for f = 2, 4, 6, . . . , v = —3, —5, —7, .... 

These constraints can be incorporated directly into the expansion by writing 
v(r,0,4>)= Yl r v b v G v {Q^)+ £ r v [a v FW(9, ft + c v F^{6^ 



(2.16) 
(2.17) 

(2.18) 



i/=2,4,6,... 
j/=-3,-5,.. 



i/=l,3,5,... 
v=-2,-4,... 



where = F, 



2G„ and fJ 2) = + (t/ + 2)G„. 
We now expand m and U + in this fashion. Requiring regularity of u at r = 
implies that a~,b~,c~ are non-zero only for v > 0. Similarly, for u + = zx + U + to 
satisfy the asymptotic condition ( 2.2 ) we must have U + — > as r — > oo and so a+, 6+ , c+ 
are non-zero only for v < 0. 

We now turn our attention to the flow on the membrane surface 8 = 1. In IHenle fcl 
|Levine (2010) it is shown that the surface flow u m (r) can be decomposed into a sum 
of incompressible shear modes w(ifj,(f>; fx) = s x V$(^>, m) where the scalar functions 
$(■0, 0; n) are eigenfunctions of the Laplacian, V 2 <I>(i/>, 4>; /x) = — fx(fx + l)$>(ip, 4>; /i). For 
the same reasons as in the bulk flow, we may take <&(ip, 0; fx) = cos (ftP^ (cos tp), whence 



w(ip,<f>;(i) = — — —V^(cos^p) cos(f)'4> ~ 2V' fl (cos ip) sin0</>. 



(2.19) 



Given a mode spectrum {/x n }, expand the velocity as u m (ij), <j>) = J2 n A n w{4>i <f>; fx n ). In- 
compressibility is satisfied for each mode by construction. Then the membrane dynamics 
plij ) become ( |Henle fc Levine||2010 1 



E 



A n [2 - Mn(Mn + l)]fln(fl n + 1)*(^> <Aj Mn) + S • V X 



2e+ 



2e; 



= 0. (2.20) 



J s=l 



There still remains the question of which modes fx to sum over. The ideal basis set 
would satisfy no-slip at ip = a in the same way as the bulk fluid basis functions. For the 
membrane basis there is not enough freedom, so we can only choose one component to 
automatically be zero. We choose to set u™(a, 4>; fx) = for all <f>, which implies 



= V' (cosa) = (1 + fx) cot aVfj, (cos a) — /icsca7 ? /J +i(cosa). 



(2.21) 



We can also take \x > by the identity V-u+i) = V^. For the case a = 7t/2, i.e. a 
perfect hemisphere (h = 0), this equation has exact solutions /i = 1, 3, 5, . . .; if a > 7t/2, 
this must be solved numerically by a method such as Newton-Raphson to obtain the 
spectrum of allowed modes. Figure [2] shows streamlines of three different recirculatory 
modes for the hemispherical case, demonstrating the increase in number of circulation 
centres as the mode number increases. The same qualitative patterns persist for h > 0. 

We now have seven sets of coefficients to solve for in the expansions: aj,6j,cj for 
u , and A n for u m . These will be determined by satisfying the simple boundary condi- 
tions (2.4) and (2.5), and the stress boundary condition (2.20). Each component of each 
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(U = 3 ;U = 5 \i = l 

Figure 2: Streamlines of recirculating modes for h = 0, projected onto the z-plane. 

boundary condition will only have <j) dependence in the form of an overall factor of sin cj> 
or cos (j), so we may factor these out and need only satisfy boundary conditions for the 
^-dependent factors. 

The complexity of the system does not lend itself easily to significant analytic progress, 
and so we choose to proceed via numerical methods to determine the coefficients. Our 
chosen solution method is the least-squares error minimisation procedure as outlined by 



Shankar (J2005J). The velocity expansions are truncated at orders N\ and N 2 for the bulk 
and membrane velocities respectively - note that use of least-squares minimisation is 
crucial here, since any attempt to solve exactly for a finite set of coefficients will yield 
poorly-determined systems. The boundary conditions are to be satisfied on < tp < a. 
Divide the interval [e,a] (e <C f) into M equal subintervals, separated by points {i^k}- 
For each boundary condition Ci to be enforced, write £\k for the error in satisfying that 
condition at ipk- Minimising E 2 — J2ik £ ik with respect to , b„ , , A n then gives a 
linear system which can be solved for the coefficients. As N ll N 2 ,M — > 00 the coeffi- 
cients quickly approach limiting values. (In what follows, wc typically use magnitudes 



N\,N2,M ~ 0(100).) Care must be taken in satisfying the stress condition (2.20), since 
the stress becomes divergent in the neighbourhood of ip — a due to the discontinuous 
geometry. To avoid this we truncate the collocation points {ipk} a small distance before 
if) = a for this condition. 

Due to the choice of origin for the bulk fluid expansion, the system quickly becomes 
numerically unstable as h increases due to divergences in r v for v > 0, r > 1 and v < 
0, r < 1 in evaluating the velocity expansions on s = 1. However, if the bulk velocity were 
expanded about the origin of the vesicle we would no longer have a basis set satisfying 
planar no-slip, and would have to enforce boundary conditions along the entire y — z = 
line as well, which would also possess highly divergent coefficients. 



3. Results and discussion 

Figure [3^, illustrates the typical flow induced on the surface of the vesicle. The closed, 
two-lobed symmetric recirculating patterns as observed by Vezy et al. ( 2007 ) are clearly 
reproduced, indicating a dynamical preference for motion rather than remaining station- 
ary. The same structure persists for non-zero values of h and for varying viscosity ratios 
r±. The recirculatory flow of the membrane induces a flow inside the vesicle appearing 
to possess entirely closed streamlines, as illustrated in figure [3f;. Similarly, the external 
flow shows regions of recirculation on either side of the membrane which also appear to 
contain closed streamlines (figure ^p) . 

The qualitative difference seen between this problem and that of an immiscible hemi- 
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Figure 3: Streamlines for r + = r_ = 0.01, h = 0. (a) Top view of membrane, (b) External 
flow, with closed streamlines in red. (c) Internal flow; brighter red indicates faster flow. 



spherical droplet (Sugiyama & Sbragaglia 2008) is due to the presence of the membrane. 
As described in f|l] simple droplet flow permits both a surface-compressible structure 
and mass exchange with the interior. Due to the membrane structure, adhesion and 
incompressibility here, no such globally circulatory flow is possible. 

The dynamical insistence on motion can be understood by the following argument. 
Consider the related problem of a free, spherical vesicle in uniform, unidirectional flow, 
for which it is known no membrane motion occurs. This problem is axisymmetric about 
the flow direction, and so the membrane experiences a tension constant along 'lines of 
latitude' perpendicular to the flow direction. Such a tension has no shearing component, 
and so the membrane is able to remain stationary. When the symmetry is broken by 
the introduction of the adhesive wall and the application of shear flow, the tension on 
the membrane is no longer uniform across 'lines of latitude'. The membrane is unable 
to support this shear, forcing it to flow. The same symmetry- breaking effects have been 
observed in vesicle membrane flows driven by electric fields ( Staykova et al. 2008 ) where 
a small field inhomogeneity develops shear stresses, causing the membrane to flow. 

It should be noted that many surface-incompressible patterns are possible whilst re- 
maining consistent with the symmetries of the problem; figure [2] shows the first three 
recirculatory modes of the membrane flow basis constructed in § |2.2| We observe that 
the flow adopts a pattern containing only two circulation centres rather than four or 
more, and pose the explanation that this is a consequence of dissipation reduction: the 
membrane is adopting the fewest number of topological 'defects' while still minimising 
dissipation within the bulk fluids. 

As the membrane viscosity is increased through the physical range, the angle i/jo 01 
the stagnation point from the vertical (see figure ^ijp) on either side shifts downwards, 
illustrated in figure [4}i. The locations seen here are close to those observed by Vezy et al. 
(2007) (where their variable 9 is equivalent to our - — ip) for small values of h (their 
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Figure 4: Dependence of the surface stagnation point angle on r± and a. (a) a = ^, 
r + = r_, 1CT 3 < r± < 10, log scale in r±. (b) r + = r_ = 1CT 2 , < h < 0.3. 




Figure 5: Linear dependence of 1/vq on l/r_ and r+, at fixed values of r + and r_, 
respectively. Top to bottom, (a) r + = 10, 5, 3.33, 2.5, 2, (b) r_ = 0.3, 0.4, 0.59, 1.11, 10. 



L/2R close to unity). We also observe that for < h < 0.3 the stagnation angle appears 
to scale linearly with the base angle a, as might be expected (figure [4J5). 

We now turn to the problem of experimentally determining the membrane viscosity 
r\ m . Due to the order of magnitude of the physical constants involved, this is a difficult 



measurement to achieve. Dimova et al. ( 1999 ) proposed a method based upon measuring 



the fluctuations of a spherical particle stuck on or penetrating through a vesicle's surface. 
While this method has been experimentally verified, and is perhaps simpler to initiate 
than encouraging vesicles to strongly adhere to substrates, it requires difficult measure- 
ments of microscopic parameters. An added benefit of the calculations presented here is 
a possible complementary method of membrane viscosity measurement based solely on 
the vesicle flow profiles, provided that the vesicle itself has not been affected by the ad- 
hesion process. Such velocities can be measured either by the defect-tracking method of 



Vezy et al. (2007) or by tracking small phase-separated patches of a second lipid species. 
Figure |5| suggests a fit of the apex velocity vq = \v m (ip = 0)| (see figure [3]d) with r± of 
the form 



"0 



Ar. 



Br 4 



Cr. 



(3.1) 



for a hemispherical membrane (h = 0). Therefore given known bulk viscosities rj±, vesicle 
radius R and measured apex velocity vq, it would be possible to solve the above relation- 
ship and so predict the membrane viscosity rj m . Using 40 x 40 regularly-spaced datapoints 
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(r + , r_) € [0.25, 10] 2 we find an excellent fit for parameter values A = 2.054, B — 4.718, 
C = 5.803, with root mean square error on the order of 10 -4 and 95% confidence intervals 
~ ±0.001. We expect such a relationship to persist for h ^ 0. 



4. Conclusions 

We have effected the theoretical solution of the problem of near-hemispherical lipid 



vesicles adhered to a flat substrate as experimentally investigated by Vezy et at. (2007 1 . 



We demonstrated the presence of symmetric recirculating regions on the membrane and 
explained the qualitative difference from the membrane-less case. Additionally, the exis- 
tence of closed streamlines both inside and outside the vesicle was shown. We then used 
knowledge of this system to suggest a new method for estimating membrane viscosities. 

Throughout this work we have assumed the membrane remains spherical. For low 
shear rates, this assumption is valid; however, behaviour at higher shear rates is unlikely 



to remain stable. Li & Pozrikidis ( 1996 ) numerically studied the deformation of adhered 



immiscible liquid drops in shear flow and found a regime where the droplet would continue 
to deform over time. It is reasonable to suggest such deformations could occur in the case 
studied here, though the more complex surface and internal flow structure may lead to 
a different parameter phase space of deformations and instability. 

Understanding the membrane flow has implications for other lipid membrane systems. 
One particular area where knowledge of the hydrodynamic behaviour of attached vesicles 
under shear flow is of importance is in the study of transient vesicle adhesion. If a vesicle 
loaded with binders comes into contact with a surface to which it may adhere, then the 
evolution of the contact area over time is dependent on the forcing experienced by the 



vesicle. Brochard-Wyart & de Gennes ( 2002 ) developed a dynamical adhesion force model 



based on the diffusion of binders on the membrane. However, if the vesicle is subjected to 
shear then free binders will be advected by the membrane flow which may significantly 
enhance or impair adhesive effects. 

Possible future work in this area includes a more complex numerical study in order to 
accommodate all sizes of spherical cap (i.e. larger h), and experimental studies using the 
proposition of ^3] to estimate membrane viscosities. It would also be interesting to study 
shear rates higher than those in Vezy et al. ( 2007[ ) to understand whether buckling and 
deformation occur. 

It is a pleasure to offer this contribution in honour of Tim Pedley's 70th birthday. We 
thank P. Khuc Trong, P. Olla and J. Dunkel for helpful discussions and S. Ganguly for 
insights at an early stage of this work. This work was supported in part by the EPSRC 
and the European Research Council Advanced Investigator Grant 247333. 
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